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BENDING AND SHEAR STRESSES DEVELOPED BY TEE INSTANTANEOUS 

ARREST OF THE ROOT OF A CANTILEVER BEAM ROTATING 

WITH CONSTANT ANGULAR VELOCITY ABOUT A 

TRANSVERSE AXIS THROUGH THE ROOT 

By Elbridge Z. Stowell, Edward B. Schwartz 
and John C. Houbolt 

SUMMARY 


A theoretical investigation was made of the behavior 
of a cantilever beam In rotational motion about a transverse 
axis through the root when the rotation of the root Is 
suddenly stooped. Equations are given for determining 
the stresses, the deflections, and the accelerations 
that occur In the beam as a result of the arrest of motion. 

The equations for bending and shear stress reveal 
that, at a given percentage of the distance from root 
to tip and at a given tip velocity, the bending stresses 
for a particular mode are Independent of the length of 
the beam and the shear stresses vary Inversely with the 
length. When examined with respect to a given angular 
velocity Instead of a given tip velocity, the equations 
reveal that the bending stress Is proportional to the 
length of the beam; whereas the shear stress is independent 
of the length. 

Sufficient experimental verification of the theory 
has previously been given In connection with another 
problem of the same type. 


INTRODUCTION 


Stresses occur in the structure of an airplane as a 
result of the shooks experienced In landing. These 
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shocks involve a rathe? Sudden ohange in the motion of 
the airplane or its component parts. 

In reference 1 a theoretical and experimental 
investigation was made of the behavior of a cantilever 
beam in transverse translational motion when the root 
of the beam was suddenly brought to rest. Equations 
were given for determining the stresses, the deflections 
and the accelerations that existed throughout the beam 
as a result of the impact. Experimental verification 
was presented for the theoretical equations. 

Shocks Involving changes in rotational motion, as 
well as shocks involving changes in translational motion 
may be experienced by the airplane. The fundamental 
theory of reference 1 has therefore been applied herein 
to an arrested cantilever beam; the only change is that 
the initial motion arrested is rotation about a trans- 
verse axis through the root. 

As in reference 1, the present paper Is based on 
the usual engineering beam theory. In this theory, the 
deflections are considered to be the result of bending 
alone and shear deflections are neglected. The theory, 
as applied to ordinary beams, gives reasonably good 
results as long as the distance between Inflection 
points is greater than a few times the depth of the beam 
When this theory for beam action Is used In vibration 
problems, such as that in the present paper, the results 
are satisfactory for those modes of vibration for which 
the nodes are not too close together. Because the 
theory of the present paper is basically the same as 
that of reference 1, the experimental verification 
given in reference 1 may reasonably be expected to 
Apply to the present equations. Additional experimental 
check was therefore considered unnecessary. 


SYMBOLS 


E modulus of elasticity 

Y weight density of material 

\ coefficient of equivalent viscous damping of 

material 
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velocity of sound In material 
acceleration of gravity 


length of beam 

moment of inertia of cross section of beam 
about neutral axis 

oross- sectional area of beam 

radius of gyration of cross’ section of 


be ami 


W) 


coordinate along beam measured from root 
distance from neutral axis of beam to any fiber 
time, zero at imoact 
operator ^ 

integers l f 2, 3* etc., designating a particular 
mode of vibration (used as subscript) 

nth positive root of 1 + cos 0 cosh 0=0 

undamped natural angular frequency of nth mode. 


radians Der second 


i-¥) , 

freqiiency of n 


damped natural angular freqiiency 

/ / 5 

radians per second l^l/l “ " 


th mode. 


When 


> 1, the "frequency" is defined 


by V = (o n| 


- 1 . 


angular velocity of beam prior to impact 
deflection of beam at station x and time t 


deflection of beam at station x and time t 
for nth mode of vibration 
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a(x,t) acceleration of beam at station x and time t 

a n ( x , t ) acceleration of beam at station x and time t 

for nth mode of vibration 

o(x,y,t) bending stress In beam at station x , distance 
from neutral axis y, and time t 

a n (x,y,t) bending stress in beam at station x, distance 
from neutral axis y, and time t for nth 
mode of vibration 

t* (X| t ) average shear stress over cross section of beam 

at station x and time t 

7 n (x,t) average shear stress over cross section of beam 
at station x and time t for nth mode of 
vibration 

A n bending- stress coefficient 

shear- stress coefficient 

C n deflection coefficient 


RESULTS AND CONCLUSIONS 
Theoretical Results 


When a cantilever beam rotating with uniform 
angular velocity 0 about the root is suddenly stopped 
at the root, theoretically an infinite number of modes 
of vibration are excited. With each successive mode, 
damping has an increasing influence upon the frequencies 
and amplitudes of vibration and, for sufficiently high 
modes, damping even changes the type of motion from 
oscillatory to nonoscillatory. In the lower modes, 
however, damping has little effect and only terms of the 
first order in damping need to be included in the 
equations. Only the equations applicable to the lower 
modes, which alone are of Importance in any practical 
case, are presented In this section. A more complete 
treatment of damping Is given In the appendix. 
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The angular frequencies (2 t r times the frequency in 
cps) are given by the equation 

0 2 

<*n = P 0 - 2 ? (1) 

IT 

where 8 n has the following values for successive modes 
of vibrations 

e x = 1 . 875104 . e 5 = 14.. 137168 

0 2 = 4.694058 e 6 = 17.278759 

e 3 = 7.854757 9 n S |(2n-l)ir, n >6 

\ = 10.995541 

These values of 8 n , which are the same as those given 
in reference 1, are characteristic of a cantilever beam. 

The energy that the beam possesses before the motion 
of the root is halted is consumed in exciting the various 
modes of vibration and is distributed among the modes as 
follows ; 


Mode , n 

Percentage of energy 

1 

97.07 

2 

247 

5 

.51 

4 

.08 

5 to pi 

.07 ■ 


This distribution of energy among the various modes of 
vibration is presented graphically in figure 1. 

Expressions for the bending stresses, shear stresses, 
deflections, and accelerations are of the same form as 
the expressions given in reference 1 for these quantities 
except that the velocity of translation v is replaced 
by the tip velocity flL and the coefficients A n , Bjj, 
and C n have new values. For the nth mode of vibration 
these expressions are as follows: 
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For stresses and deflections. 


= A n ® 


sin <0 n t 




■■sir' 


sin a> n t 


» . » _ fiL It 

V*.*) = c n — 7 


sin u) t 
n 


For accelerations, when damping is sufficiently small, 

a n (x»t) = -u> n 2 w n (x,t) 

The variations of the dimensionless coefficients A n , 
B n , and C n with position along the beam x/L are 

given for the first three inodes (n = 1, 2, and 3) in 
figures 2 to Ij.. The highest absolute values of A n 
and Bjj, and hence the highest stresses, occur at the 
root of the beam. These values, for the first three 
mode s , are 


Mode , n 

A n at root 

B n at root 

1 

1.1376 

I.566 

2 

. iBiS 

.8679 

3 

.06I4.O 

.5088 


The foregoing values of A n and B n at the root are 
presented graphically in figure 5- 

VThen the effects of damping are neglected, the 
maximum values with respect to time of o n (x,y,t) and 
“ n (x,t) associated with the nth mode of vibration are 


o n (x,y) = A n ^ % 
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and 


T„(x) = bJ£ 


n' 


c L 


( 7 ) 


Equations (2) to (5) for stresses, deflections, 
and accelerations give the values associated with the 
nth mode of vibration. Since all modes of vibration 
occur simultaneously, the net results are the super- 
position of the effects of all modes. Hiis super- 
position gives the following equations; For bending 
stress. 


o(x,y,t) ■=— Z. eIais 

o p \ 


2E 


\(D . 


BintOjt + Age 


2E 


sin ogt + 


For average shear stress. 


.. Ol p 
tT x, t) * — “■ E e 
o L 



sin (i^t + B2© 


Xcc 2 2 

2E 


sinWgt + 


( 6 ) 


(9) 


For deflection. 


-ai£ 

o p 



Xw^ 

2E 


sin + Cge 


Xd>g 2 
2E 


Bin Wgt + 


(10) 


For acceleration, when damping is sufficiently small. 


a(x, t) 


QL 
o p 



Xd > p 2 

t 1 

Bind) t + C^0g 2 e sin d)t + 



01) 


'Che equations for bending stress reveal that, at 
a given percentage of the distance from root to tip 
and at a given tip velocity, the bending stress for a 
particular mode is Independent of the length of the beam. 
The equations for shear stress reveal that, for a given 
tip velocity, the shear stresses at any station vary 
inversely with the length of the beam. When examined 
with respect to a given angular velocity instead of a 
given tip velocity, the equations for bending and shear 
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stresses reveal that the bending stress is proportional 
to the length of the beam; whereas the shear stress Is 
Independent of the length of the beam. 


Influence of Type of Shock 

The frequencies and deflected shapes of the modes 
of vibration are characteristic of the cantilever and 
not of the Initial shock. The type of shock, however, 
determines the relative amplitudes of the various modes. 
The amplitude of each mode. In turn, determines the 
energy and stresses developed by that mode. In fig- 
ure 6, arrested translation (reference 1) and arrested 
rotation (present paper) are compared on the basis of 
the distribution of deflections, energies, and stresses 
among the first three modes. In the case of arrested 
translation, about 6l percent of the total energy Is 
in the first mode; the first mode predominates as 
regards tip deflection and root bending stress, while 
the largest shear stress occurs In the second mode. In 
the case of arrested rotation, about 97 percent of the 
total energy Is In the first mode, and the first mode 
predominates not only with respect to tip deflection 
and root bending stress but also with respect to 
maximum shear stress. 

If the motion that Is arrested is a combination 
of translation and rotation about the root, the deflec- 
tions and stresses can be found by superposing the 
results for the translation and the rotation separately. 

If the cantilever is not perpendicular to the direction 
of root translation at the moment of Impact, the 
translational velocity used should be the component of the 
root velooity perpendicular to the cantilever at the 
moment of Impact. The component of root velocity parallel 
to the cantilever causes longitudinal vibrations but 
such vibrations are not considered herein. 

The equations given in this report and In refer- 
ence 1 apply also to the cases In which the root of 
a cantilever beam at rest Is suddenly given a constant 
angular velocity or a constant linear velocity in a 
direction perpendicular to the length of the beam. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX 

THEORETICAL D’ERIVATIOTT 1 ' 


General analysis .- Consider a beam of uniform cross 
section in equilibrium. If a portion of the beam la 
suddenly disturbed, as by a shock, the beam Is set Into 
damped bending oscillations. The equation of motion 
for these bending oscillations Is given by the dif- 
ferential equation (reference 1) 



+ \p2 


d^w , Y = 0 
dx^- ‘ dt g £t^ 


(Al) 


By using 


p gn 

c = — , equation (Al) can 
Y 

dKi + \ &5 W + d 2 w 

E dx^ dt p^c^ dt^ 


be written 


= 0 


(A 2 ) 


This partial differential equation Is reduced to an 
ordinary differential equation of the fourth order by 
d 

writing p = — -j thus 



(A 3 ) 


The general solution of equation (A 3 ) Is 

w = P cosh 6 7 + Q sinh 6~ + R cos 87 + S sin 87 ( aJU- ) 

L Ii L L 


where 


8 = 



ip 


P°\/ 1 + P g 


The coefficients P, Q, R, and S are to be determined 
from the boundary conditions. The case under consideration 
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is that of a cantilever beam rotating about its root 
with constant angular velooity ft and having the 
rotation suddenly stopped by clamping of the root in 
oosition as the beam coincides with the x-axis. The 
boundary conditions for this case are 




0 


The angular velocity of the root as given by the 
boundary condition 0 - QJ? (discontinuous function) is 
represented graphically in figure 7. By the procedure 
adopted in reference 1, the solution will be obtained 
for the boundary condition 



= -M 


and the constant angular velocity O will be added to 
the resulting angular velocity. 

With the application of the boundary conditions to 
equation (aIj.), the operational solution for the angular 
velocity induced by -Cll is found to be 



2 ( 


- . ’ft (i 

1 + oo8h 6 oos 8) 


0)(o 


+ 008 6 oosh 0)(oo8h0j- + 008 0 


i) 


^oob 0jr “ oosh 0j^ 


+ sin 0 sinh 0 /oos 0r- - cosh 0 
+ (oosh 0 sin 0 - sinh 0 oos 


0)^*inh0£- + sin0~^ 


Interpretation of this operational expression by use of the Heaviside expansion 
theorem and addition of the constant angular velocity ft gives for the total 
angular velocity . 


d/w 

dtVfcXi 


00 

}= ft - fti +ft ^ p, { 0 nf)‘ 


\w n 2 


cos <*> n f t - 



sin o) n »t fA5) 


nth positive root of 1 + cos 0 cosh 0=0 (all roots - 
1 namely, ±0„ and ±10 n - have been considered in the 

_ interpretation) 

0 n 2 

a = pc — p- undamped natural angular frequency of nth mode. 

It- radians/seo 

I x2 *»n 2 

a » = Wqv /I - damped natural angular frequency of nth mode,. 

V ljJS 2 radians/sec 

/ v sin 0 n sinh ©nfcos 0^ - oosh 0 n £) + (oosh 0 n sin 0 n - sinh 6 n oos 0 n ) (slnh0 n “ + «in 

>'(0 nr) -2 S ^ 7 — ' 1 r — 1 

' / 0 n yjosh 0 n sin 0 n - sinh 0 n oos 8 n j 


Integration of equation (A5) with respect to x gives for velocity 
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Integration of equation (a 6 ) with respect to time with the condition (w) t _ Q 
gives for the deflection (when t ^ 0 ) 


w(x, t ) = QL 


^2- F(0 7 ) ^-t 

/ — * e sin w n 'ti 

n=T V 




where 
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The contribution of the nth mode to the deflection 


w n (x,t) = ~ “ c n 
c P 



sin “ n »ti (a8) 


When > 1, equation (A8) may be put in the form 

2E 


/ . » _ fill It „ 

w n (x,t) - — — c n 


\0) n 2 


/\2 <dn 2 


sinh w n * 1 1 (A9) 


where now 


J n* “ u ny 1^2 


The form indicated by equation (a8), in which - — < 1, 

is characteristic of the lower modes and represents 
damped oscillatory motion. The form indicated by 

\o> n 

equation (A9), in which — > 1 (damping greater than 

cE 

critical), is characteristic of the higher modes and 
represents subsidence motion. 


The comule te behavior of the cantilever may be 
determined from equations (a 6) for velocity and (A7) 
for deflection. The quantities of interest are the 
bending stress, the shear stress, and to some extent 
the accelerations. When damping is present, the 
equations representing the contribution of the nth 
mode to these quantities may be given in the two forms 
indicated by equations (a 8) and (A9). In subsequent 
equations, however, only the form indicated by (a 8) is 
given, because this form is characteristic of the modes 
that are of practical importance. 


Bending stresses.- The bending stresses a(x,y,t) at any fiber located a 
distance y from the neutral axis are 


QO 


\o>. 


o(x,y,t) = Ey 
where 


= J = = = 

4x2 ° ,/ris 

V / V ktf. 


n 


2E gin a> n f t^ 


An " 2 


T * 

^oosh 8 n aln 8 n - alnh 8 n coa 9^^oosh9 n £ + °°* e n jr) - 9 ° Jlnh 9 n( aljh9 nf ♦ alllfl nf) 


©n^^osh Qq^ sin 0 U - sinh 0 n cos B n ^ 


Tbe contribution to bending stress of the nth mode is 


= E ~ ^ - 1 


\o n 2 


r l - 


e 22 sin &> ’ t J? 


kst 


Average shear stresses.- The average shear stress on the cross section is 


T( x ,t) = Ep 2 = E — £ 
3 c L 


dx" 



Xw n 2 


Bn 


s= e 2E sin w 'ti 


l^ 2 


n 


1 
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where 


Bn -2 


■ 2i 


cosh 6 n Bin‘0 n - sinh 


0 fl oog 0^ |slnh 6^— - sin 0 - sin 9 n sinh 8 n ^QBh 6^2.+ oos % ^ 


0 n ^ccsh 0 n sin 9n - sinh 0Q oos 3 n ^ 
The contribution to shear stress of the nth mode is 


>oo. 


n 


T n U,t) = E ^ ^ 

c L n 


/l-^b 

lj.E 2 


2E sin co n *t 1 


Accelerations .- Prom equation (a 6), with the aid of the relation 

pF(t)i = F(0)pi + F* ( t )J? 

the acceleration anywhere on the beam is found to be 
a(x, 



e 


^n 2 

2E 



sin oo n »t + 


1 - 


— COS 0)_ 1 1 \1 

XAon 2 


2E 2 


H 

Vi 
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With the aid of the orthogonal properties of the function 
to show that 

p (Mh=° 

0=1 

The linear accelerations are therefore always finite. 

The contribution to acceleration of the nth mode is 

e “ I sin 0 ) 


Comparison with the expression for w n (x,t) from equation (a 8) shows that the 
acceleration for each mode is out of phase with the deflection. When damping 
is sufficiently small, however, the relation between the acceleration and the 
deflection reduces to the well-known result for undamped vibration 

an(x,t) = -w n 2 w n (x,t) 


n 


»t + 


Xun* 

E 


2 - cos V t 


\ 0) 


1 - 


n 


2E 


a n (x, t ) =- 


fiL i£ 


n w n 


1 - 


1 - 


x2(0 n 2 

2E 2 

" A ? 


Iie 


,2 


'(M)* 


it is possible 


On 
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Fig. 5 
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Figure 5.- Values of bending -stress coefficient An 
and shear- stress coefficient B n at root, = O . 
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(a) Deflection at tip. 




(c) Maximum bending stresses at root. 



(d) Maximum shear stresses at root. 


I Translation 


fl Rotation 


Figure 6. - Comparison of principal results of reference I 
(cantilever in translation) and the present paper 
(cantilever in rotation) . 
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Figure 7. - Graphical representation of the 
discontinuous function CL-CII. 





